Theory of Relcixation Dynamics in Glass-Forming Hydrogen-Bonded Liquids 



H.G.E.Hcntschcl^''^ and Itamar Procaccia^ 
^Dept. of Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 
^Dept of Physics, Emory University, Atlanta, Georgia. 

We address the relaxation dynamics in hydrogen-bonded super-cooled liquids near the glass tran- 
sition, measured via Broad-Band Dielectric Spectroscopy (BDS). We propose a theory based on 
decomposing the relaxation of the macroscopic dipole moment into contributions from hydrogen 
bonded clusters of s molecules, with Smin < s < Smax- The existence of Smax is due to dynamical 
arrest and its value may depend on the cooling protocol and on the aging time. The existence of 
Smax is translated into a sum-rule on the concentrations of clusters of size s. We construct the 
statistical mechanics of the super-cooled liquid subject to this sum-rule as a constraint, to estimate 
the temperature-dependent density of clusters of size s. With a theoretical estimate of the relax- 
ation time of each cluster we provide predictions for the real and imaginary part of the frequency 
dependent dielectric response. The predicted spectra and their temperature dependence are in ac- 
cord with measurements, explaining a host of phenomenological fits like the Vogel-Fulcher fit and 
the stretched exponential fit. Using glycerol as a particular example we demonstrate quantitative 
correspondence between theory and experiments. The theory also demonstrates that the a peak 
and the "excess wing" stem from the same physics in this material. The theory also shows that in 
other hydrogen-bonded glass formers the "excess wing" can develop into a /3 peak, depending on the 
molecular material parameters (predominantly the surface energy of the clusters). We thus argue 
that a and (3 peaks can stem from the same physics. Finally we address the BDS in constrained 
geometries (pores) and explain why recent experiments on glycerol did not show a deviation from 
bulk spectra. 



I. INTRODUCTION 
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In this paper we propose a theory, based on a sim- 
ple physical model at the mesoscale, that can account 
quantitatively for most of the observed features in the re- 
markable relaxation dynamics of glass-forming hydrogen- 
bonded liquids [11, i, i, i, i, i, Q . The paradigmatic ex- 
ample of such systems is dry glycerol and glycerol-water 
mixtures, but other alcohols and alcohol mixtures show 
similar properties. The relaxation dynamics of such sys- 
tems in the vicinity of the glass transition exhibit an ex- 
tremely wide range of frequencies, spanning sometimes 
16 orders of magnitude or more in the frequency domain. 
Despite the fact that good measurements have been avail- 
able for more than half a century and after consid- 
erable experimental and theoretical effort, including ap- 
parent "first principles" theories like Mode-Coupling 3, 
there is no accepted derivation of the observed charac- 
teristics of the relaxation dynamics. Until now most dis- 
cussions of experimental data are limited to fitting phe- 
nomenological exp ressions that are not the result of a 
proper theory [lO|. The present paper attempts to close 
this gap. 

An excellent technique to probe this enormous range of 
frequencies is broad-band dielectric spectroscopy (BDS) 
that covers the frequency band (10~^- 10^^ Hz) and 
temperature range (100-600K). The interpretation of the 
measurements in BDS is facilitated by the fact that the 
dielectric theory for BDS is well established and inde- 
pendent of the nature of the relaxational mechanisms 
involved in any particular material and at any partic- 
ular temperature and pressure. We can express the fre- 
quency dependent dielectric constant e(w) in terms of the 
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Here is the high frequency dielectric constant due to 
fast rotational processes above lO^^Hz and atomic polar- 
ization; eo is the static dielectric constant; and adc is the 
dc conductivity of the medium. The response function 
(j){t) is given by the normalized correlation function of 
the macroscopic dipole moment M(t) 
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A typical BDS absorption spectrum read from the 
imaginary part of e(tj) with dry glycerol at temperature 
T — 196K is shown in Fig. [TJ One observes what has 
been termed in the literature "the main relaxation pro- 
cess" or a peak at fmax ~ 10~^Hz, and then tapering off 
initially with a typical power law form uj~^ which will 
appear as a straight line on a log- log plot, but is gen- 
erally followed by a typical wing, which in Fig. [T] has a 
gentle curvature at frequencies above / > lO^Hz. This 
part of the spectrum had been termed in the literature 
" the Excess Wing" ; This term does not reflect a deep re- 
alization of the existence of an excess loss, but rather the 
name stems simply from the inability of the phenomeno- 
logical Davidson-Cole fit formula [ll| to agree with the 
experimental spectrum at this frequency range. While 
not in glycerol, in other hydrogen-bonded glass formers 
this wing can also exhibit a clear shoulder or even a dis- 
tinct second (3 peak ^ . The literature does not agree on 
whether the a relaxation and the 'excess wing' (or the 
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FIG. 1: A typical BDS absorption spectrum read from the 
imaginary part of e((jj) for dry glycerol at temperature T — 
196. The dashed line is the phenomenological Davidson-Cole 
formula, which fails to describe the so called "excess wing" . 
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FIG. 2: Experimental BDS absorption spectra of dry glycerol 
at a range of temperatures 202-292 K, where frequencies were 
rescaled to the maximum of the a peak and amplitudes were 
rescaled to the amplitude ot he same peak. Note that in 
comparison to Fig. [T]here the dc contribution is shown. The 
data collapse indicates strongly that the a peak, the "excess 
wing" stem and the low frequency conductivity contribution 
are due to the same physical mechanism. 



/3 peak when it exists) stem from the same physics or 
not 01 . A very strong experimental argument in favor of 
the very same physics was presented in Q where spec- 
tra taken at different temperatures were rescaled by the 
position and amplitude of the a peak, see Fig. [2j The 
excellent data collapse is a strong indication that the a 
process and the "excess wing" stem from the same phys- 
ical mechanism. The theory presented below will corrob- 



orate this finding and will clarify the physical origin of 
the a and P peaks when both exist. 

Of course, the raw spectra exhibit a strong tempera- 
ture dependence. Using spectra measured in the range of 
temperatures 190 — 2A0K in dry glycerol, the positions 
of the maxima of the "main process" have been fitted to 
the Vogel-Fulchcr formula by defining T„iax = ^/'^'^fmax 
and then writing 



T„iax{T) = exp [DTu/ {T -Ty) 
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The three parameters were fitted to the data with the re- 
sults Inr^, = -35.9, D = 22.7 and T„ = 125K. Note that 
this fit implies an attempt time Ty w 10~^^ sec which 
is a very short time indeed. We will argue below that 
the Vogel-Fulcher formula has no deep meaning, and the 
parameters involved are of limited physical interest. Al- 
though the spectra derived from our theory below will 
be shown to agree with Eq. ([3]) in the range of measure- 
ment used, we argue the formula is nothing but a data fit 
that can be used only in a finite range of temperatures. 
The prediction of the theory described below is that that 
there is no true divergence of T,nax at any temperature 
T > 0. 

The task of the theorist is then to derive, preferably 
from first principles, the form of the frequency dependent 
dielectric response e{uj) for this system, to explain the ob- 
served spectra (both real and imaginary parts) . In partic- 
ular a theory should indicate whether the "main process" 
and the "excess wing" stem from the same physics or does 
one need to invoke more than one relaxation mechanism. 
We will see that the former is the case here. Then one 
needs to explain the observed fits to the 'Vogel-Fulcher 
formula. While the present authors do not ascribe much 
physical significance to the parameters appearing in the 
formula, nevertheless the experimental fits should be ac- 
counted for, and hopefully the origin of the parameters 
identified. Finally one needs to explain the data collapse 
shown in Fig. [51 We believe that the present paper 
achieves these tasks. 

In order to understand the broadband dielectric spec- 
troscopy of glycerol (and in the future water-glyccrol mix- 
tures and other alcohols), we will treat glycerol in the 
relevant range of temperatures as a heterogeneous fluid 
on macroscopic timescales. That is, that while on very 
long timescales the liquid phase must be homogeneous, 
there exist localized mesoscale domains in the fluid that 
have macroscopic lifetimes. Indeed, inhomogeneities that 
appear to survive for 10^ seconds contributing to the di- 
electric response in the Fourier domain at frequencies as 
low as 10~^ Hz in some cases. We shall develop the 
theory on the basis that these inhomogeneities are a dis- 
tribution of clusters having the structure of an incipient 
strongly hydrogen-bonded phase and a surrounding bath 
of more mobile and less dense "liquid-like" glycerol phase. 
Together therefore a snapshot of this phase will present 
disordered "glass-like" characteristics. Our task is two- 
fold: First, to express the dielectric response in terms of 
this cluster distribution; and second, to find the distribu- 
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tion of clusters. We shall do this using mesoscale thermo- 
dynamic arguments. Once we have these distributions we 
describe the resulting BDS, and we shall see the appear- 
ance of both a dominant "a" peak at low frequencies and, 
depending on molecular parameters, an 'excess wing' or a 
secondary shoulder at higher frequencies. We should 
note that in spirit our approach combines ideas based 
on the molecular dynamics observations by Geiger and 
Stanley [l^ on the appearance of low-density patches 
in hydrogen-bonded water (though in the present con- 
text we are actually considering higher density patches in 
hydrogen-bonded glycerol), together with ideas of Cham- 
berlin [l^ and Kivelson et al. [13, [IBl on treating such 
heterogeneties as clusters. The similarity to the latter 
theories extends however only up to Eq. ([6]) below, af- 
ter which our theory diverges from theirs. The crucial 
differences from previous theories will be pointed out be- 
low as we go along. We trust that the reader will find 
the present approach superior in simplicity and in the 
quality of the predictions. 

In Sect, [n] we present the dielectric theory for a 
medium in which there exists clusters whose dipole mo- 
ments are responsible for the dielectric response. We 
estimate the relaxation time of such clusters, and most 
crucially, their concentrations. In computing the con- 
centrations we deviate most strongly from all previous 
theories. In Sect. Illll we specialize the theory to the case 
of glycerol; it is interesting to see how, by using a mod- 
est input of experimental data we get naturally a host 
of results in very good correspondence with experiments. 
We explain the apparent Vogel-Fulcher and stretched ex- 
ponential fits, but also point out that the range of va- 
lidity of these fits is finite, whereas the present theory is 
more widely applicable. Sec. IIVI opens up the discussion 
for other materials, and examines the changes expected 
in the spectra when molecular parameters change. We 
show how the relative prominence and a and /3 peaks 
depend on such parameters. In Sect. |V]we discuss BDS 
in constrained geometries and explain why glycerol in 
small pores exhibited the same spectra as bulk glycerol; 
the explanation is that the smallest pores were just big 
enough to contain the largest clusters predicted by the 
theory. Reducing the pores a bit further should result in 
a dramatic change in the spectra. Section [VII presents a 
summary and a discussion. 



II. DIELECTRIC THEORY FOR CLUSTERS 

A. The Physical Model and the Fundamental 
Assumption 

We propose that the long-time relaxation in glycerol 
arises due to the existence of clusters of s molecules 
which are hydrogen bonded to make the cluster a rec- 
ognizable entity. We denote the number of clusters of 
s molecules as Ng. In terms of these clusters we can 
write the dipole moment of the system, M{t), as the 



sum M{t) = J2a T^a{t) + Miiguid{t) whcrc nia{t) is the 
dipole moment of cluster a and the fast relaxing fluid 
contribution to the dipole moment Muquidit) will not be 
studied in detail as its relaxation spectrum is expected 
to be significant at high frequencies. We expect that in 
the glassy state relaxation phenomena are relatively rare 
events, allowing us to assume that different relaxation 
events are statistically uncorrelated. This is the funda- 
mental assumption of the model, i.e. that the relaxation 
process in each cluster is statistically independent of the 
other clusters. In fact, by a "cluster" we mean a bunch 
of molecules that are highly correlated; if there are two 
adjacent clusters that are highly correlated, they should 
be considered as one cluster. With this in mind we can 
then write 
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We can rewrite Eq. [5] in terms of an intensive number 
density Us = Ns/M , where M is the total number of 
molecules in the system: 
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To each cluster we associate its longest typical relaxation 
time T, which we identify with the cluster lifetime. Thus 
(f){t) assumes the form 

"■s(ms - m,) exp (— i/r,) , , 

Mt) = 1^ , ' — ^\ ' ' . 6 

2^^n,(m, - m,) 

To proceed further we need to estimate n,, (m, -m,) and 
T,. We start with the static dipole correlations. 



B. The static dipole correlations 

The dipole moment of a cluster of size s can be ex- 
pressed in terms of individual glycerol molecule dipoles 
dj as m, = X]i=i ^i- F'or a set of non-interacting dipoles 
(m, - m,) = s(P , where d = \di\. We do not expect 
the dipoles in the clusters to be totally random, but 
rather to exhibit strong short-range correlations due 
to dipole-dipole interactions. This short range order was 
taken into account explicitly by Kirkwood who intro- 
duced the so-called Kirkwood factor g: (m, -m,) = gs(P, 
where g > 1- Obviously this constant is difficult to cal- 
culate for a given material, but luckily it appears in both 
the numerator and denominator of (f)(t) and consequently 
Eq. ^ reduces to the expression 
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If we insert Eq. [7] into Eq. [T] and split the dielectric con- 
stant into its real and imaginary parts e(w) = 5Re(tj) -I- 
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iQe{uj) we find 
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Thus we see that we have expressed both the broadband 
dielectric constant and loss in terms of the cluster size 
distribution and the cluster lifetimes. 



C. Relaxation times Ts 

The relaxation time of the clusters, mainly due to their 
rotational relaxation, will be determined by the free en- 
ergy barrier that involves breaking the hydrogen bonds 
with the surrounding liquid. We argue that these are 
given by Arrhenius forms where the energy barrier scales 
with the surface area of the cluster as the clusters at- 
tempts to break the cage of mobile 'liquid-like' molecules 
in which it is confined. 



Tac^j>ijls^/^/kBT) 
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Here the attempt time tq is of the order of 10^^^ seconds; 
while the energy for breaking a typical bond ji k, a can be 
expected to scale with the surface energy per molecule. 
We will estimate the numerical values of these parameters 
below in the context of the theory for glycerol. 

While Eq. ^ appears rather innocent and perfectly 
reasonable, its consequences are manifold, lying at the 
very basis of our approach. These consequences are ex- 
plained in the next subsection. 



D. Mesoscale Thermodynamic Theory for the 
Cluster Distribution ns{T,p). 

The most interesting part of the theory is the esti- 
mate of the cluster size distribution ns{T,p). In this 
section we provide a thermodynamic theory to estimate 
this crucial quantity. Obviously, any such approach has 
to face the fact that the glassy state is thermodynam- 
ically metastable, with the ground state always being 
the crystalline solid, i.e. a single "cluster" with s ~ 00. 
The actual nature of this metastable state depends on the 
preparation protocol. As is well known, different rate of 
cooling for example can end up in glassy states where the 
properties of the metastable state are different. More- 
over, left alone, such a state will "age" @, changing ever 
so intermittently and ever so slowly towards the ground 
state. What we are really talking about is dynamical 
arrest. 

To see the role of dynamics, note that if the relaxation 
time as expressed in Eq. ^ is correct we can estimate 



the size of the largest clusters as Smax just from looking 
at the longest relaxation times observed experimentally. 
We thus write 



(10) 



where tmax is the longest time scale observed in the 
broadband spectrum. For example, for glycerol tmax ~ 
10'' seconds and p. ^ 2kBT; then Smax ~ 100. Of 
course, this number can depend on the cooling proto- 
col, and when the system ages we expect Smax to be a 
slowly (and probably intermittently) increasing function 
of time. Nevertheless wc will see below that in order of 
magnitude this estimate is fully justified by the present 
theory. But more importantly, we can now turn the ar- 
gument around, and say that if the cluster distribution is 
stationary, then the relaxation time of the cluster should 
be of the same order as its creation time. So let us es- 
timate, if it takes 10^ seconds to create a cluster with 
Smax = 100, how long will it take to create a cluster of 
size 2smax- Using Eq. (jTU]) with tq « 10""'^^ yields the 
estimated time to create a cluster of 200 molecules to be 
longer than 10^^ seconds! This is the crux of the matter. 
The dynamical arrest can be interpreted here as having 
a distribution of clusters in which Smax is sharply de- 
termined by the aging time of our glassy state. In the 
language of phase transitions we can say that the phase 
transition had occurred, but the coarsening of the clus- 
ter distribution is inhibited by dynamical arrest. In any 
reasonable theory of coarsening, like, say, the Lifshitz- 
Slyozov theory [13|, the time scale is determined by the 
diffusion coefficient, which in our super-cooled liquid is 
many orders of magnitude smaller than in regular liq- 
uids undergoing a phase transition. Thus at any time we 
have a maximal cluster, and the smaller clusters which 
form much quicker have had ample time to equilibrate 
according to the regular laws of thermodynamics. 

Thus, In developing a thermodynamic theory, we need 
to invoke a constraint on the size of the maximal clus- 
ter. In other words, we treat the heterogeneities in glyc- 
erol in terms of a mesoscale distribution of clusters with 
a range of size Smin < s < s„iax in a more mobile 
bath of 'liquid-like' molecules. This approach of applying 
to the glass-forming liquid statistical mechanics subject 
to constraints follows up o ur p revious analysis of other 
glass forming liquids, cf. |l8l. flQl. [20| . Note that this 
approach is radically different from the proposition of 
[ij, [ist where a mysterious "strain energy" term was in- 
voked to explain why larger clusters are not present in 
the theory. 

The intensive number of such clusters of size s are 
ns{p,T) = Ns{p,T)/AI, and the probability that an in- 
dividual glycerol molecule belongs to a cluster of size s 
is Cs = nsS. We denote below the number of molecules 
belonging to the clusters and to the liquid phase by Mc 
and Mi respectively, with Mc + Me = AI. To describe the 
mesoscale thermodynamics we assume that in a cluster 
at temperature T and pressure p we have an intensive 
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contribution to the chemical potential 

lJLc=Uc+pVc-TSc, (11) 

where the subscript c stands for "cluster" and Uc, Wc, 
Sc are respectively the internal energy /molecule, the vol- 
ume/molecule, and the entropy/molecule in the clusters. 
Similarly in the 'liquid-like' phase we can write for the 
intensive contribution to the chemical potential 

= ue+ pve - Tse, (12) 

where the subscript £ stands for liquid and Ui, ve, se 
are respectively the internal energy /molecule, the vol- 
ume/molecule, and the entropy/molecule in the mobile 
'liquid-like' phase. We expect that Uc < ug because of 
the the hydrogen bonds in the clusters are less distorted; 
Sc < sg because of the greater number of rotational de- 
grees of freedom in the mobile phase; and Vc < vg be- 
cause of the higher density of more solid-like glycerol 
compared to the liquid. All these intensive thermody- 
namic variables are in principle functions of temperature 
T and pressure p. Apart from these extensive contri- 
butions to the Gibbs free energy Mcfic + Mgfif where 
Mc = J2s = Mj2s ~ ^^J2s "^s' need to con- 
sider two other crucial contributions to the total free en- 
ergy. 

First, there is a surface energy contribution U surf ace 
due to the interface between each cluster and its sur- 
rounding bath. This scales as the surface area of each 
cluster and therefore we can write 

C/sWace ='T^7V,s2/3, (13) 

s 

where a is the surface energy /molecule. Obviously, the 
value of a will turn out to be crucial in determining the 
distribution of cluster sizes. To estimate its size note 
that a could be as large as a hydrogen bond energy, but 
in fact we expect it to be smaller because the mobile 
phase is also partially hydrogen bonded. Thus we expect 
a ~ Ui ~ Uc which is some fraction of a hydrogen bond 
energy. 

Second, there will be an entropic contribution due to 
all possible spatial configurations of clusters in the glassy 
liquid. This is an important contribution to the free en- 
ergy, favoring small clusters and being responsible 
for the /3 peak when it exists. In taking the en- 
tropy explicitly into account we deviate once more from 
all previous theories; we cannot see why the entropy was 
not considered before. We can estimate this entropy by 
sequentially placing all clusters in the available volume 
starting with the largest [I^ . Thus the volume available 
to the largest clusters is y = McVc + M^vi. The volume 
available to the next largest clusters i&V — SmaxNg^^^Vc 
etc. Proceeding in this manner the entropy of mixing is 
given by the expression 

Smix = -ks'^Mslxs Inxs -I- (1 - Xs) ln(l - x^)] (14) 



where fc^ is Boltzmann's constant and Ms = {V/vc — 
Tlis'>s ^s's')/ s is the number of boxes available to clus- 
ters of size s while Xs = Ns/Ms are the fraction that are 
occupied. Thus 

Xs = ^F^s —j 1 ^ • (15) 

where ci is the fraction of molecules in the 'liquid-like' 
phase, and of course the sum of the concentrations of 
molecules in the "condensed' and 'liquid-like' phases obey 

Cc + Q = 1.0. (16) 

At this point we introduce v as the average volume per 
glycerol molecule in our system, where of course 

V = CcVc + ClVi . (17) 

We can now solve Eqs. ((TB)) and (fT7|) for Cc and cg in 
terms of these 3 volumes per molecule and write the cru- 
cial sum-rule: 

> Cs=Cc= . (18) 

^ Vc - Vi 

S — Sjjiin 

This is the crucial constraint on the thermodynamic the- 
ory ; we will impose this constraint on our solution by 
choosing the largest cluster s = Smax such that Eq. (fT5|) 
is obeyed. 

The Gibbs free energy for our system is given by the 
expression 

G + PV -TS Mc^ic + Mefie + a ^ NsS^^^ 

S 

+kBT"^Ns[xs Inxs + (1 - Xs) ln(l - x,)]. (19) 

S 

and the chemical potential of a cluster of size 5 is given 
by 

Us = sfi = dG/dNs = fics + (Ts^/^ 
+kBT{\nijs + Y.^s/s')\n{\-'i>s')} (20) 

where 

ijs = X ' (21) 

Vg Vc — V 

A = — — . 

Vc Vc - Vi 

Note that < A < 1. Rewriting Eq. [20]for the chemical 
potential per molecule in cluster 

H = fic+<7S-'^^ + kBT{{l/s) ln^, + ^(l/s')ln(l-i^s')} 

s'>s 

(22) 
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The chemical potential in the liquid phase must be the 
same as that calculated for the clusters and given by 
Eq. [22l and therefore 

fi^fie + kBT{vi/vc){{l/s) In + ln(l - Vs)}- 

(23) 

We can use Eq. (|23p to define a new chemical potential 
/i' = /i—(/.i—yU£)/(wf/wc)—Mc which together with Eq. ([22]) 
then obeys 

^' = as-'/"" + kT{{l/s)ln[4's/{l-0-i^s)] 

- 5](l/s')ln(l-^.0} ■ (24) 

s' <Cs 

In this form it is easy to solve these equations numerically 
in a sequential fashion starting with s ~ Smin- Using 
Eq. ([2T|) we see that the concentration of the smallest 
clusters Cg^-^ obeys 

= ^-C.™' + kBT{l/s„-,„,) ln[c,„„„ /A] . (25) 

To estimate Smm we note that in 3-dimensions a 'clus- 
ter' must have at least 2-3 molecules in each direction to 
be identifiable as a cluster. Thus the minimal value of s 
is about 20. We do not know ji' a priori depending as 
it does on the details of the energies, structure and en- 
tropies in the clusters and 'liquid-like' phases, so instead 
we choose the reasonable values for fi' ~ ksT. We are 
now in a position to estimate Cs„i„ from Eq. 1251 Then 
all the values of Cg for s > Smin follow from Eqs. ((M]) 
and (HI]) in order until Eq.([THl) is fuUfilled. We will see 
below that for glycerol the largest s predicted by the the- 
ory is about 100. 

III. PREDICTIONS FOR BROADBAND 
DIELECTRIC SPECTROSCOPY: THE CASE OF 
GLYCEROL 

Examining our expressions for the distribution of clus- 
ters given by Eqs. we note that in order to 
apply the theory to a particular substance, we miss the 
value of the surface energy a, say for glycerol. Since 
this quantity appears inside exponential forms we need a 
rather accurate value for comparison with experimental 
data. One approach would be to fit our dielectric spectra 
at one temperature and then use the fitted parameters 
to predict its form at different temperatures. Here we 
have taken a slightly different approach, using available 
glycerol data. The are two types of data; the first con- 
cerns the average volumes per glycerol molecule in the 
liquid, solid and glassy phases. For Vi we employed data 
of liquid glycerol at room temperature giving vi = 1214"^. 
For Vc and v we took data from [2lj which give for the 
the cluster phase Vc = llOA^ and for the glassy phase 
V = 11 3 A'^ at atmospheric pressure. These numbers were 
used to estimate the constraint given by Eq. The 



second type of data we have at our disposal concerns the 
size of dynamic heterogeneities as a function of temper- 
ature. We have taken data from for the average size 
of such heterogeneities Sav = J2s ^'^s/ J2s '^s- The data 
used were Sav{T = 210i4:) « 50, Sav{T = 205 A') « 58, 
Sav{T = 200^:) « 67, Sav{T = 195K) w 73. These 
two pieces of information can be accommodated in our 
theory using cr/fes = 320iir or a/ksT w 1.6. We then 
predict the maximal cluster sizes at these temperatures 
to be s„^aAT = 210K) « 63, s„,,(T = 205K) « 72, 
SmaooiT = 200/i) « 81, s,naociT = 195K) « 87. Using 
this data we then studied the theoretical predictions for 
glycerol. 

As noted above, the distribution of clusters between 
Smin and Smax IS strougly dependent on surface energy. 
Using a/ks = 320i^, we plot the cluster size distribution 
ris and the density of molecules as a function of cluster 
size s for dry glycerol at T = 210K, 205K, 200A", 195A^ 
Small clusters are favoured at higher temperatures, when 
the entropy of mixing dominates the distribution (see 
Fig. [3]). But as the temperature is reduced larger clusters 
are favoured as energy dominates the distribution. In 
general the distribution is bimodal favouring small and 
large clusters at the expense of intermediate sizes (see 
Fig. . We note that the sharp cut of the distribution 
at s = Smax is a bit artificial; in reality one can expect 
a sharply decaying tail at s values slightly larger than 
Smax- We do not expect such minor details to influence 
the main results presented below. 

We are now in a position to calculate both (j>{t) from 
Eq. ([7]) and the real and imaginary parts of the dielec- 
tric function from Eq. ([8|). This quantity had been fit- 
ted phenomenologically to a stretched exponential form, 
(j){t) ^ exp —{t/r)^'^, a form referred to in the literature 
as the "Kohlrausch- Williams- Watts (KWW) relaxation 
fimction". (cf. for example [1^). To see whether this 
form is justified by the present theory we plot our com- 
puted function in the appropriate coordinates, cf. Fig. 21 
Indeed, a stretched exponential with 0.4 < fSx < 0.66 
gives an acceptable fit over a broad range of timescales 
tmin 10~^ <^ t <^ tmax ~ 10° secs, with deviations 
at both shorter and longer times. We note that Pk is 
temperature dependent and stress that the stretched ex- 
ponential form has a limited value in the sense that it is 
just an acceptable fit in a limited range for a very dif- 
ferent function. Nevertheless the numerical value of Pk 
and its variability with the experimental conditions are 
both confirmed by experiments in glycerol, see for exam- 
ple gm. 

Of greater interest is the dielectric spectrum and loss 
function (see Fig. [5|) which correspond to the distribu- 
tions shown in Fig. [3] with T = 210K, 205K, 200K, 195K 
respectively). The real part of the dielectric constant 
^e{iu) start to decline from eo to its asymptotic value 
of Eoo at around uj « Wmax but takes several frequency 
decades to achieve asymptotia. We refer the reader to 
[1, [2^ and references therein to note that both the qual- 
itative form of 5Re(cL')and its quantitative details are in 
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^ cluster size distribution for glycerol at T=200K 



Stretched exponential plot of i|)(t) versus log .(t) 
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FIG. 3: Upper panel; Cluster size distribution versus s for 
temperatures T = 210K, 205K, 200K, 195K. Lower panel: 
Concentration of molecules Cs as a function of cluster size s 
for temperaturesT = 210K, 205K, 200K, 195K. 



close correspondence with experiments in glycerol (as 
well as in glycerol-rich water mixtures). 

The shape of 3e(u;) is controlled by the shape of the 
cluster size distribution. There is always a clear a peak 
at low frequencies uJmax defined by d^e{u!rnax) / di^max = 
associated with the largest clusters. Very roughly 
i^max ~ ^/tmax- But wc Can also sce at all tempera- 
tures a clear 'excess wing' at higher frequencies. As the 
temperature is lowered the 'excess wing' becomes slightly 
less pronounced. 

The next interesting question is whether we can jus- 
tify theoretically the Vogel-Fulcher plots. To answer this 
question wc plot the computed values of LUmax in a Vogcl- 
Fulcher form as in Eq. cf. Fig [HI In this plot we 
display, in purpose, all the temperature range includ- 
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FIG. 4: Stretched exponential plot of (f>{t) showing and ap- 
proximate straight line over about six orders of magnitude 
in time for temperatures T = 210K (leftmost), 205K, 200K, 
195K (rightmost). 



Real part of dielectric constant for glyceroi divided by its static vaiue 
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FIG. 5: Upper panel: Real part of the dielectric spec- 
trum e'{uj)/eo versus log]^Q(tj) for temperatures T = 
2WK, 205K, 200K, 195K. Lower panel: A log- log plot of the 
dielectric loss logjQ(e"(ci;)) versus logj^g(aj) for temperatures 
T = 210K (rightmost), 205K, 200K and 195K (leftmost) 
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Volger-Fulcher plot for glycerol 
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FIG. 6; ujmax for the dielectric loss for temperatures T — 
210/«:, 2057^, 2007^, 195/r plotted in a Vogel-Fulcher form. 
The plot underlines the irrelevance of T„ as a physical pa- 
rameter, since it is so far removed from the region of linear 
fit. The theory presented here actually does not predict a 
true divergence of the relaxation time, it only becomes super- 
exponentially large. 



ing to stress the absurdity of such a fit. Neverthe- 
less a straight line can be fitted through the computed 
points. Using the same (unphysical) attempt time of 
lnr„ = —35.9 as in the experimental fits on glycerol the 
straight line best fit gives = 129K and a fragility 
D = 20.0. This should be compared with = 125K 
and D = 22.7 which are the numbers reported experi- 
mentally. Our conclusion is that the theory explains both 
the stretched exponential fits to the relaxation function 
and the Vogel-Fulcher fit, but both are fundamentally 
meaningless, and should be replaced by a theory of the 
type proposed here. 
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IV. THE ROLE OF SURFACE ENERGY: 
CLUSTER DISTRIBUTIONS AND DIELECTRIC 
SPECTRA 

To understand the experimental results for glycerol we 
chose the crucial parameter a from structural data. In 
this section we ask a different question — what is the 
qualitative influence of the surface energy a on cluster 
size distributions and dielectric responses. In particular 
we are interested in the possibility of generating a dis- 
tinct P peak by changing only one molecular parmater, 
which is a. To study this question we fix Smax = 100 and 
T = 200K and study the effect of changing the surface 
energy. We choose four values a/ksT = 1.0, 1.25, 1.5 
and 1.75. As we increase a the average cluster size in- 
creases: Savicr/ksT = 1.0) w 48, Savifr/ksT = 1.25) « 



FIG. 7: Upper panel: Cluster size distribution Us versus s for 
a/ksT = 1.0, a/ksT = 1.25, a/ksT = 1.5, cr/fcsT = 1.75. 
Lower panel: Concentration of molecules Cs as a function 
of cluster size s for temperatures a/ksT = 1.0,a/fcsT = 
1.25, a/ksT = 1.5, a/ksT = 1.75. 



76, Savi<y/kBT = 1.5) w 87, Say{a/kBT = 1.75) « 91. 
The complete distribution can be seen in Fig. [71 As tr 
increases large clusters are favoured over small clusters; 
(7 has a crucial qualitative infiuence on the bi-modality 
of the distribution. 

The change in surface energy must also infiuence the 
dynamics. We asserted that the lifetimes Ts of clusters 
of size s are given by Arrhenius forms where the energy 
barrier scales with the surface area of the cluster, as the 
clusters attempts to break the cage of mobile 'liquid-like' 
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Stretched exponential plot of (t)(t) versus log (t) 

10r 




l°9io« 

FIG. 8; Stretched exponential plot of to show scaling 
regime in time for a/ksT — 1.0, a /knT — 1.25, a /ksT — 
1.5, a/feflT = 1.75. 



molecules in which it is confined. The energy for breaking 
a typical bond can be expected to scale with the surface 
energy a. In Fig. [8] we have plotted the behaviour of 
(j){t) in a stretched exponential form to show scaling in 
time for a/ksT = 1.0, 1.25, 1.5 and 1.75. Note that 
as a increases the stretched exponential regime increases 
due to a more pronounced a peak while at low surface 
energies the high frequency (3 peak destroys this kind of 
scaling. 

Of greater interest are the dielectric spectra and loss 
functions (see Fig. ^ which correspond to the distribu- 
tions shown in Fig. [7] with a/ksT = 1.0, 1.25, 1.5 and 
1.75 respectively). The shape of 5fte(w) corresponds to a 
decline from eg to eoo starting at tOmax, but one interest- 
ing point is that the (3 peak can lead to a shoulder here 
also. The shape of 3e(tt') is again controlled by the cluster 
size distribution. There is always a clear a peak at low 
frequencies. But for lower values of the surface energy 
small clusters are encouraged resulting in a prominent 
(3 peak at high frequencies. As the surface energy in- 
creases the /3 peak becomes less pronounced turning first 
into a shoulder at intermediate frequencies (see Fig. [9]), 
and finally into the anomalous scaling observed in many 
experimental data sets (see Fig. [9]). These interesting 
qualitative findings will be turned to quantitative com- 
parisons with experiments in different materials in a later 
publication. 
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FIG. 9: Upper panel: Real part of the dielectric spec- 
trum e'(ti^)/eo versus logjQ(cj) for a/ksT = 1.0 (rightmost), 
a/ksT = 1.25, a/kaT = 1.5, a/ksT = 1.75 (leftmost). 
Lower panel: A log-log plot of the dielectric loss logjQ(e"(a')) 
versus logjQ(ctj) for the same values of a/ksT as in the upper 
panel. 

V. BDS IN CONFINED GEOMETRIES 

Finally, we discuss the BDS of glycerol in confined ge- 
ometries. We refer here in particular to the experimental 
studies reported in [2^ , in which it was shown that glyc- 
erol could be confined in pores whose diameter d can be 
as small as d — 2.5 nm in size without seeing any ap- 



10 



Frequency peak for glycerol at T=200K 
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FIG. 10: Shift in expected position of ujmax as a function of 
confining pore diameter. Note that no effect appears until a 
critical pore diameter after which the shift is dramatic. 

preciable effect in the position of the a maximum or its 
amplitude. Since we expect that the maximum size clus- 
ter should be limited by the size of pores, it appears sur- 
prising that there is no confinement effect on the spectra. 
If our arguments are correct as presented in this paper, 
there must be a point at which the dielectric loss will be 
strongly effected by confinement. To estimate where this 
should happen we equate the volume of the pore to the 
volume of the largest cluster, i.e. 

TTcf/G « S,naxVc ■ (26) 

Using our estimates of Vc and referring to the temper- 
atures employed in [23 |. we find that confinement ef- 
fects are expected to appear when the diameter satis- 
fies d < 2.5 nm. Thus the experiment had just missed 
the confinement effect by a hair. In Fig, [10] we present 
our own prediction as to how the expected position of 
ijJmax depends on the confining pore diameter. Note that 
there is no effect down to a critical pore diameter below 
which the shift in cumax is dramatic. Here we are plotting 
shift as a function of pore diameter at fixed temperature. 
Wc also note that as a function of temperature clusters 
tend to grow or shrink, and therefore confinement effects 
should be more dramatic at lower temperatures where 
the unconfined clusters are expected to be larger. 

VI. DISCUSSION 

We have shown how a relatively simple model of un- 
correlated cluster relaxations can account, without much 
parameter-fitting, for the qualitative and even the quan- 
titative aspects of the observed BDS in hydrogen-bonded 



Scaled plot of dielectric loss for glycerol 
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FIG. 11: Data collapse with the theoretical spectra: arescaled 
version of the plots in Fig. [5] to which the dc conductivity 
contribution to the loss has been added. 

liquids. The theory was demonstrated for glycerol, but 
obviously with a mere change of molecular parameters it 
should apply to a broad range of other hydrogen-bonded 
liquids. It was demonstrated explicitly that the a and 
(3 regions of the dielectric spectra can in principle stem 
from the very same physics, and their relative amplitudes 
is determined by the relative population of small or large 
clusters. The relative population is determined by en- 
tropic effects which were explicitly taken into account in 
our thermodynamic theory, and by the surface energy per 
molecule in addition to the temperature. Applications to 
other materials will be presented elsewhere. 

One particularly interesting aspect of the glycerol mea- 
surements, i.e. the data collapse shown in Fig. [21 was 
not used yet to challenge the theory. Note that the ex- 
perimental data collapse includes the clearly identified 
dc branch below frequencies f / fmax ~ 10""^. We present 
our theoretical spectra in a similar manner, rescaling the 
frequency uj and the amplitude e"{uj) by the frequency 
and amplitude of the a peak. The result of the exercise 
is shown in Fig. 1111 without any attempt to re-fit any of 
the material parameters. Note the approximate data col- 
lapse, with the excess wings failing to collapse perfectly. 
It is clear that by choosing an appropriate value of a 
we could achieve a better data collapse, but wc do not 
try to do so here, deferring detailed quantitative fits to 
a future publication where theory and experiments will 
be compared in full detail. We stress at this point that 
the dc branch was added here by hand, and we do not 
have a clear understanding why it is involved in the data 
collapse together with the a peas and the 'excess wing'. 
The existence of such scaling at low frequencies suggests, 
however, that the dc conductivity can be written in the 
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form Udc ^ ^max^max- Why this is so requires further 
investigation. 
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